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ABSTRACT 

The ubiquity of M dwarf stars combined with their low masses and luminosities 
make them prime targets in the search for nearby, habitable exoplanets. We investi¬ 
gate the effects of starspot-induced radial velocity (RV) jitter on detection and char¬ 
acterization of planets orbiting M dwarfs. We create surface spot configurations with 
both random spot coverage and active regions. Synthetic stellar spectra are calculated 
from a given spot map, and RV measurements are obtained using cross-correlation 
technique. We add the RV signal of an orbiting planet to these jitter measurements, 
and reduce the data to “measure” the planetary parameters. We investigate the de¬ 
tectability of planets around M dwarfs of different activity levels, and the recovery of 
input planetary parameters. When studying the recovery of the planetary period we 
note that while our original orbital radius places the planet inside the HZ of its star, 
even at a filling factor of 2% a few of our measurements fall outside the “conserva¬ 
tive Habitable Zone”. Higher spot filling factors result in more and higher deviations. 
Our investigations suggest that caution should be used when characterizing planets 
discovered with the RV method around stars that are (or are potentially) active. 

Key words: stars: activity, atmospheres, planetary systems, rotation, spots, low- 
mass - techniques: radial velocities 


1 INTRODUCTION 

M dwarfs are of great interest in the search for nearby hab¬ 
itable planets for many reasons. They have been hailed by 
some as the holy grail of stars for potential life-permitting- 
planet hosts, but have been condemned by others as un¬ 
suitable targets for planet searches. Due to their ubiquity 
alone -M dwarfs make up ~70 percent of stars in the Milky 


Wav (Henrv et al. 

1994j; Chabrier 2003j; |Reid et all [2004: 

Covev et all 20081: 

Bochanski et al. 2010j), combined with 


the fact that observations indicate the occurrence rate of 
super-Earths and Earth-sized planets in creases with decreas¬ 
ing stellar mass llHoward et al.l boid l , the discovery of a 
nearby habitable planet seems most li kely to be found o rbit¬ 
ing an M dwarf. A recent study bv lDressing fc Charbonneaul 
d2013l) found a lower limit on the occurrence rate of Earth- 
size planets in the habitable zones (HZ) of cool stars of 0.04 
planets per star, and also determined that the nearest planet 
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in the ha bitable zone should be just 5 pc away, although 
iMann et al.l lj2013l f contests that this estimate might be too 
high based on their revised values of R *, L„, M* for a sam¬ 
ple of Kepler stars, which resulted in some of the formerly 
HZ planets no longer orbiting within the HZ. 

The HZ is most loosely defined as the area around a 
star in which liquid water could exist on the planet’s sur¬ 
face. A ccording to the model developed bv lKopparapu et al.l 
Uni), the inner and outer HZ limits for our Solar System 
are 0. 99 AU and 1.70 AU, respectively. iKopparapu et all 
(1201311 point out that their model does not include the ra¬ 
diative effects of clouds, so the boundaries could extend far¬ 
ther. M dwarfs are much less massive and less luminous than 
solar-type stars, which means the HZ is significantly closer to 
the star (for example, for an M dwarf of mass M ~ 0.5Mg, 
and R ~ 0.5Rq, the HZ is between ~0.1 and ~0.2 AU). This 
makes small planets orbiting in the HZ around M dwarfs eas¬ 
ier to detect than their counterparts orbiting more massive 
stars for a few reasons: (1) Transits are more likely to be 
detected due to the the HZ being significantly closer to the 
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star lKastinelll993h . which results in a larger range of in¬ 
clination angles that will yield a visible transit from Earth. 
A HZ planet orbiting a 0.5 M@, 0.5 Rq M dwarf has a ~3 
times higher likelihood of transiting compared to a planet 
orbiting a Solar-mass star in the star’s HZ. (2) A larger 
planet-to-star ratio results in larger transit depths in the 
light curve, since the transit depth is simply proportional to 
the area of the star that is being eclipsed by the transiting 
planet. An Earth-mass planet in the HZ of a R = 0.5 Rq M 
dwarf would have ~4 times the transit depth compared to 
a similar planet orbiting a Sun-like star. (3) The close-in 
habitable zone around M dwarfs also makes radial velocity 
(RV) reflex motion larger (and thus in principle easier to 
measure). Planets of a given mass would be easier to detect 
around M dwarfs than around larger stars because of the 
reduced mass ratio of the system. An Earth-mass planet or¬ 
biting in the HZ of a 0.5Mq M dwarf would have a Doppler 
shift 3 to 4 times stronger than the same planet orbiting in 
the Sun’s HZ. 

However, a HZ near to the star can cause a host of other 
potential problems. For instance, a planet orbiting close to 
the star can result in tidal heating from the planet itself 
(I Barnes et al.l l2013l l, as the gravity from the star squeezes 
the planet and creates inner friction. Additionally, many 
M dwarfs are high l y active, with lon g activity lifetimes 
(|Hawlev et al.l l l 99fll : IWest, et al.l 120081 4 . Magnetic activity 
from the star can potentially render even a planet well within 
the habitable zone inhospitable to life. Flares can irradiate 
the surface of the planet, or dramatically affect the atmo¬ 
sphere, possibly removing it completely (see for example Se¬ 
gura et al. 2010). 

Stellar activity can also inhibit the detection of orbiting 
planets. Dark spots on the surface of the star—a result of 
magnetic activity in the photosphere—can create noise in 
light curves, and radial velocity “jitter”: a term which refers 
to the noise-like structure introduced into the measured ra¬ 
dial velocity curve in a systematic or unsystematic way, for 
example when the wavelength of a stellar line from the in¬ 
tegrated stellar surface light is shifted, during the stellar 
rotation, due to the existence of inhomogeneities on the sur¬ 
face intensity caused by sunspots, plages, or other phenom¬ 
ena moving semi-permanently with the rotating surface or 
appearing and disappearing during the rotation. In first ap¬ 
proximation one could assume that this phenomenon could 
only look like a planet with an orbital period identical to the 
stellar rotation period. However, such jitter in the data can 
cause several different kinds of misinterpretations of the real 
signal, including false-positive planet identifications, mask¬ 
ing the planetary signal, or leading to an erroneous estimate 
of the planetary parameters. The jitter from stellar activity 
has the largest probability of dominating the Fourier peri- 
odigram when the planetary signal is weak relative to the 
stellar activity signal, i.e. for small planets and/or stars with 
large jitter-creating activity, and/or when the obtained data 
are too sparsely sampled in time. 

Spots that persist over timescales significantly longer 
than the stellar rotation period can mimic a planetary sig¬ 
nature, resul ting i n false-positiv e “detections” of planetary 
companions. iQueloz et all d200lf) examined the photometri¬ 
cally variable and magnetically active star HD 166435, which 
had displayed low-amplitude RV variations with a period 
of just under 4 days, implying a close-in planetary com¬ 


panion. However, closer investigation revealed that the RV 
signal w as not coheren t berond durations of approximately 
30 days. iQueloz et alJ (j 200lf ) concluded that the RV varia¬ 
tions could be explained by line-profile changes due to dark 
photospheric spots on HD 166435, and were not due to a 
gravitational interaction with an orbiting planet. 

There are various methods used in attempts to distin¬ 
guish between spot- and gravitationally-induced RV vari¬ 
ations. Unfortunately, no consistently reliable method has 
been found. Bisector analysis is one common method of at¬ 
tempting to classify the origin of observed RV variability. 
Spots will distort the absorption line profiles from the star 
and result in a change in the asymmetry of the lines. These 
distortions—and the resulting asymmetries—change as the 
star rotates and the spot crosses from one stellar limb to the 
other. The distortions can be quantified with the line bisec¬ 
tor span—the difference in bisector value from two different 
locations in the line profile. The size of the apparent RV shift 
resulting from the spots is expected to be correlated with the 
size of the line bisector span. Unfortunately, bisector analysis 
is not always usefu l in distinguishing the cause of RV vari¬ 
ability. IPrato et al.l d2008T ) demonstrated a case where young 
stars with RV variations that likely resulted from spots did 
not have significant correlations between these RV varia¬ 
tions and the bisector spans. Also, if usini is lower that the 
resolution of the spectrograph, variations i n RV and bisector 
spans will not correlate fsee lDesort et al1l2007l f. 

The effects of starspots on RV m easurements have been 
studied by e.g. iMakarov et al.l 1 2009 1. who derived the am¬ 
plitude of astrometric, photometric, and RV perturbations 
caused by a single spot, to compare the sensitivity of astro- 
metric method to the Doppler technique. They concluded 
that if the ultimate limit of exoplanet detection is defined 
by intrinsic astrophysical perturbations, the astrometric 
meth od is more sensitive t han the Doppler technique in this 
limit. iBoisse et al.l (j2011al l performed simultaneous model¬ 
ing of stellar activity and planetary parameters in an at¬ 
tempt to remove the RV jitter. They succeeded in removing 
up to 90% of the spot-induced jitter from their test models, 
but required well-sampled photometry and accurate charac¬ 
terization of the host star in order to distinguish th e jitter 
from the planet-induced RV signal. iDumusaue et all d2014l l 
developed the software tool Spot Oscillation And Planet 
(SOAP), which simulates the effect of stell ar spots and 
plages on radial velocimetry and photometry. iBarnes et ali 
(1201 lT l studied how spot-induced jitter affects planet studies 
of M dwarfs by using randomly-distributed sp ot maps to ob ¬ 
tain d etection thresholds for orbiting planets. IReiners et al.l 
d20ldl investigated spot-induced jitter as a function of obser¬ 
vational wavelength as part of a larger study of detection of 
pla nets around very low-mass stars using RV measurements, 
and lMarchwinski et ah ( 2014 ) used the Sun as a proxy to in¬ 
vestigate RV jitter as a function of wavelength by estimating 
the RV jitter induced by stellar activit y on the Sun us i ng so - 
lar spectral irradiance measurements. iReiners et al.l (2010) 
noted that the jitter depends strongly on the properties of 
the spot, including spot temperature and appearance. 

All of these studies are affected by a limited knowl¬ 
edge of the behavior, amount, and patterns of spots on 
stellar surfaces, especially those of M dwarfs. Often for M 
dwarfs, solar models are extrapolated to more active stars, 
and random/uniform spot coverage is assumed, which is not 
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necessarily correct. Doppler Imaging allows “mapping” of 
starspots on the stellar surface, but its application to M 
dwarfs is limited due to the need for fast rotators combined 
with bright magnitudes, the latter of which limits the num¬ 
ber of M dwarfs that can be observed. New techniques such 
as Eclipse Mapping through exami nation of starspot cross- 
ing events in photometric data fsee ISanchis-Oieda fc Winnl 
Ecmj) have revealed more information about spot param¬ 
eters on M dwarfs stars and make possible exploration of 
these parameters beyond the limits of Doppler Imaging. 

In this work we combined information from Doppler 
Imaging, Eclipse Mapping, light curve monitoring, TiO ob¬ 
servations, and dynamo models to paint a more realistic 
picture of M dwarf spots, including spot coverage, spot-to- 
photosphere temperature contrast ratios, and active regions 
on the stellar surface, which we discuss in Section 2. Section 
3 explains the methods we use to generate line profiles and 
extract the stellar RV signals. We first probed the variation 
in jitter with changes in model parameters such as obser¬ 
vational wavelength, stellar temperature, and spot config¬ 
urations. We then added a simulated planet to the model, 
and briefly explored the planet detection limits imposed by 
stellar jitter. In Section 4 we investigated the uncertainty 
contributed by stellar jitter to derived planetary system pa¬ 
rameters. Finally, our resulting limits on exoplanet detec¬ 
tion and characterization due to stellar spot-induced jitter 
are discussed in Section 5, as well as implications for planet 
searches, and future work. 


2 M DWARF SPOTS 


Unfortunately, despite the ubiquity of small stars, little is 
still known about M dwarf spots (due, in part, to their low 
luminosities). We attempted to model correct spot behavior 
based on current known constraints and spot information 
from the literature. 

One manner in which stellar spots are parameterized 
is by quoting the spot filling factor, the percentage of the 
visible surface area of the star that is covered by dark spots. 
Spot filling factors can be measured in various ways, from 
direct observations in the case of Sunspots (see ISolankil 
19 9ll), to o b servat ions of absorption bands of TiO (e.g. 
O’Neal et al. 2004ll . and to some extent, Doppler Imaging 
(see iRic S)7 although this technique will not resolve 
small spots and thus will underestimate filling factors. Spot 
filling factors can also be indirectly inferred by creating 
model stars with various spot sizes and filling factors and 
comparing the generated light curves to observational data. 
However, simply knowing the spot filling factor is not suffi¬ 
cient, since different configurations of spots will yield signif¬ 
icantly different amounts of RV jitter. An extreme example 
of this would be to imagine two stars, each with a near 50% 
filling factor, where on one star the spots were uniformly 
distributed across the stellar surface, and on the other star 
all the spots were grouped on one hemisphere. These two 
scenarios would result in vastly different observational sig¬ 
natures. A less extreme version of this scenario is illustrated 
in Figure [I] where three simulated stellar surfaces are shown. 
Each has the same spot filling factor, but a different size 
distribution of spots, resulting in a differing level of homo¬ 
geneity of the spot coverage for each star. The derived RV 
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Figure 1. Top: Three stellar surfaces with the same photosphere- 
to-spot temperature contrast (AT = 500) and spot filling factor 
(30%), but different spot size distributions. Each was generated 
by randomly ch oosing spots from a lognormal distribution of spot 
radii (following IBogdan et al1ll988l f. combined with a power law 
distribution to increase the number of small spots. The mean 
value, //., for the lognormal distribution is the only factor that 
changes, and corresponds to (a), fi r = 0.250 deg, (b), fi r = 0.700 
deg, (a), fi r = 1.500 deg. Bottom: RV jitter generated by the 
three different spot maps, at 20 observational phases spanning 
the entire longitude of the stellar surface. 


jitter from each surface shows how even a small difference 
in spot size can result in a significant (in this case, a factor 
of almost 3) change in spot-induced jitter. 

To correctly model M dwarf spots, we made use of con¬ 
straints on spot size, temperature, location, and behavior 
from various observational methods such as Doppler Imag¬ 
ing, Zeeman Doppler Imaging, investigation of light curve 
modulation, and eclipse mapping, as well as modeling of M 
dwarf dynamo behavior. Evidence for both random, uniform 
spottedness as well as spots concentrated in active latitudes 
and/or longitudes has been observed. 


2.1 Spot size and distribution 

Doppler Imaging (DI) operates on the same principle as 
spot-induced jitter (dark spots on the stellar surface result 
in an anomaly in the observed spectrum), by “inverting” 
the spectrum to extract the surface distribution of spots 
that created it. This is one of the few ways to obtain a 
“map” of a stellar surface, although the solution is degen¬ 
erate, so results should be taken with caution. Due to their 
low-luminosities, fewer than ten K and M dwarfs have been 
successfully imaged using the DI technique. Of these, the 
results are mixed: 

Some low-mass stars have shown evidence of high- 
latitude (or polar), long-lived spots or spot groups, that 
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appear at some epochs and are not present at others. For 
example, AB Dor, a rapidly-rotating, ZAMS K dwarf dis¬ 
playe d this type of spot behavior in 1993-1994 (| Jeffers et al.1 
l2007lf . but not in 1989 dKuerster et al .1 [l 9941) . and BO Mic 
displayed similar high-latitude spots in 200 2 dBarnesll2005l l 
but not in 1998 dBarnes fc Cameronl l200lh . This behavior 
could be indicative of stellar cycles on these stars. 

The rapidly-rotating M dwarfs HK Aqr (Ml, P=0.4 
days) and EY Dra (M2, P=0.5 days) show spots at 
low latitudes or with no strong latitude dependence 
l|Barnes fe Cameronl l200ll ). though based on photometric 
observ ations EY Dra_has shown indication of active longi¬ 
tudes (IVida et al.ll20ldh where spots preferentially appear. 

Further evidence for non-uniform spot distribution is 
provided by the Eclipse Mapping technique, where proper¬ 
ties of starspots can also be inferred through observation of 
anomalies in high signal-to- noise, high-caden c e light curves 
of planetary transits (see Pont et al.l 120071 : IWolter et al.l 
120091 : ISanchis-Oieda fe Winnl 201 ill . This has some advan¬ 
tage over Doppler Imaging since the stars do not need to 
be rapidly rotating and it is possible to obtain informa¬ 
tion about less luminous M dwarfs than would be possible 
through Doppler imaging. In eclipse mapping, the eclips¬ 
ing body “scans” the face of its companion and passes over 
dark spots, and the resulting changes in light level can be 
used to recons truct the surf a ce bri g htness distribution of the 
eclipsed face. Wolter et a.1.1 (l2009l 'l. ISilva-Valio et alJ (l2010h 
and lSilva-Valio fe Lanzal ( 201 il l found numerous Spot Cross¬ 
ing Events (SCEs) in the light curve of CoRoT-2, a rapidly- 
rotating (P=4.5 d) G7V star. They found that an average of 
five SCEs occurred per transit, and that the transited region 
had a spot filling factor of 10 20 percent. 

Eclipse mapping is still difficult for M dwarfs because 
it requires a transiting object to scan the face of the star. In 
the K epler catalog, there are ap proximately 4000 M dwarfs 
([Dressing fc Charbonneaull2013l) . Of these, so far only one 
has been identified that has both a transiting planet and a 
corresponding anomaly in the light curve indicative of an 
SCE: KOI 1152. The dearth of M dwarf SCEs could im¬ 
ply either low contrast ratios between spot temperatures 
and photospheric temperatures, or stars that have high fill¬ 
ing factors of more or less homogeneously distributed spots 
(or simply that the probability of a planet crossing a spot¬ 
ted region is low). KOI 1152 indicates a large spot size 
of approximately 20 degrees in diameter, and a low tem¬ 
perature contrast ratio of about 100 degrees difference be¬ 
tween the spot temperature and the photospheric temper- 
ature of the star (Sanc his-Ojeda, private communication). 
ISanchis-Oieda fe Winnl (l201lh also identified spot-crossing 
anomalies on the K4V star HAT-P-11 which preferentially 
occurred at two specific phases of the transit. They inter¬ 
preted these as intersection points between the transit chord 
and active latitudes of the host star. Based on these data, it 
would make sense to model both homogeneously distributed, 
high filling factor spot configurations as well as large spots 
and/or spot concentrations in specific active regions. 

Analysis of the broad-band light curves of magnetically 
active stars seems to yield indirect evidence of a more ran¬ 
dom, uniform distribution of spots. The amount of rotational 
modulation in the stellar light curve depends on the degree 
of axisymmetry of the starspot distribution. Thus, observa¬ 
tion of a lack of rotational modulation in the light curve 


could be indicative of either a spotless surface or a very uni¬ 
form spot coverage that does not significa ntly alter the lumi¬ 
nosity of the visible surface with rotation. I Jackson fe Jeffriesl 
d20l3) examined magnetically active low-mass members of 
NGC 2516. While some displayed small amplitude modula¬ 
tions with rotation, they noted that almost half of the mon¬ 
itored stars had no detectable light-curve modulation at all, 
despite showing similar chromospheric emission levels (in¬ 
dicative of magnetic activity), having a similar distribution 
of equatorial rotation velocities, and showing no significant 
difference in distribution of positions on a color-magnitude 
diagram. They suggested it should be possible to have high 
spot filling factors and still display small light-curve ampli¬ 
tudes if the stars were covered by many small, randomly 
placed spots with angular diameters of ~ 2°. 

Ijackson fc Jeffries] d2013l J further examined this pro¬ 
posal through simulations of low-mass magnetically ac¬ 
tive stars with randomly orientated stellar spin axes and 
cool starspots of a characteristic scalelength randomly dis¬ 
tributed across the stellar surfaces. They compared the mod¬ 
els with their observations of NGC 2516 M dwarfs and found 
that the best-fitting starspot angular scalelength is about 3.5 
degrees, although this is dependent on the the assumed spot 
temperature ratio and filling factor (i.e. the model cannot 
be constrained by the light curve data alone). 

ISavanov fc Dmitrienkd (120121 '). however, found evidence 
of two active longitudes on the fully-convective M dwarf LHS 
6351 using light curve monitoring, similar to the I Vida et al.l 
d2010h active longitudes. 

Zeeman Doppler Imaging (ZDI) seems to further sup¬ 
port the idea of a more uniform spot distribution. ZDI al¬ 
lows mapping of the magnetic topography of the star by 
detecting spectral line splitting due to the Zeeman effect. 
Like regular Doppler Imaging, the result is degenerate, and 
the reconstruction of the magnetic geometry yields only the 
simplest solution that fits the observed data. There is evi¬ 
dence from ZDI of more axisymmetric magnetic field topol¬ 
ogy inMuhyj^convectiyestars (approximately M4 and later, 
see iMorin et alJl2008l. l2010h . 

Dynamo simulations of M dwarf stars have been carried 
out, though the resulting consequenc es on spot forma tion 
(and location) are often not mentioned. iBrownind d2008h ran 
three-dimensional nonlinear magnetohydrodynamic simula¬ 
tions of the interiors of fully convective M dwarfs to ascer¬ 
tain if a magnetic dynamo could be produced without the 
shear caused by the interface between c onvective and ra - 
diative zones as in larger stars. Although Br ownin 3 (|2008|J 
concluded that magnetic fields can be generated in fully- 
convective stars, and that the fields would have a strong 
axisymmetric component, they stopped short of speculating 
on the resulting spot formation on the stellar surface. 

iGranzer et al.l (l2000l l performed simulations of the dy¬ 
namics of magnetic flux tubes in low-mass ZAMS stars (with 
masses between 0.4 Mq and 1.7M@), and calculated the spot 
emergence latitudes at different rotation rates. They found 
that for a 0.4 Mq star, spots form at least 20 deg from the 
equator even at the slowest rotation rates. At higher rota¬ 
tion speeds, they found spots concentrated between ~50 to 
70 degrees in latitude. 
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2.2 Spot temperatures 

Spot radius and temperature are highly correlated when fit¬ 
ting light-curve (LC) and eclipse-mapping (EM) data. Thus, 
a knowledge of spot temperatures would be ideal, although 
difficult to obtain. Spot temperatures can also be obtained 
using molecular band head modeling (MB), line depth ratios 
(LDR), and DI. DI seems to indicate that cooler stars show a 
lower temperature difference between photosp h eric temper¬ 
ature and spot temperature (IStrassmeierll2011ff . iBerdvuginal 
(l2005ll showed an average temperature difference for M 
dwarfs of just 200 K, while for larger stars AT values were 
obs erved up to 2000 K . 

iGrav fc Johansonl (Il99lh demonstrated the LDR tech¬ 
nique of using two nearby spectral lines of different tem¬ 
perature sensitivity to obtain spot tempera t ures. This 
method was emplo yed by e.g. ICatalano et, al.1 (120021 ') and 
iFrasca et al.l (l2005h who obtained a range of 400-1000 K 
for AT for a number of giant stars. 

Due to their low temperatures, cool dwarfs exhibit 
many TiO bands which can be u sed to derive spot tem¬ 
peratures, using the MB technique. lO’Neal et alJ (2004) ob¬ 
tained relative spot temperatures for five G and K dwarfs 
using the TiO band heads at 705 nm and 886 nm and found 
AT values in the range of 1000-2000 K. 

iRice et ahl d2010l ') obtained two separate sets of Doppler 
images of the WTTS V410 Tauri star, from atomic lines and 
from the TiO 705.5-nm lines respectively, and found that 
the TiO-based temperatures appeared to be cooler, with an 
average temperature difference of 105 K between the two 
techniques. This could partly be a result of the difficulty 
in determining the correct continuum in the forest of TiO 
lines, however they concluded that atomic-line DI is prone to 
underestimate spot temperatures on a level of about 100 K 
for a 4500 K photospheric temperature. 

Figure [2] shows spot temperature contrasts from the lit¬ 
erature, including both giant stars (Tablc[lJ and dwarfs (Ta¬ 
ble [5J. Only twelve data points (for eight individual stars) 
exist for M dwarfs. We attempted to fit a line to these points 
but the resulting line suffered from small numbers, and had a 
slope opposite to the linear and polynomial fits that include 
more points, likely indicating that more data are needed be¬ 
fore such a fit will realistically reflect the change in spot 
contrast with photospheric temperature of M dwarf stars. 
Linear fit lines are also of limited use for predicting temper¬ 
ature contrasts at the low photospheric temperatures of M 
dwarfs due to the fact that the lines reach zero in the T = 
2400 - 2600 K range. Since the data are so sparse and the 
error bars so large, it seems the most we can currently con¬ 
clude about M dwarf spot temperature is that AT seems to 
be in the range of 200 - 600 K, with a couple stars showing 
higher contrasts of up to ~800K. 


3 SIMULATING RADIAL VELOCITY DATA 

We use these observations and simulations of spot behav¬ 
ior as constraints on our modeled M dwarf spot patterns. 
We create cases with random, “uniform” spot coverage 
(i.e. no preferential latitude or longitude spot emergence) 
as well as cas e s wit h act ive latitude region s following the 
iGranzer et al. | (l20Q(]ll and lVida et al.l d 200 oT ) models. 



Figure 2. Spot temperature contrasts from published values (see 
Tables [l] and [2j. Filled points indicate dwarfs and open points are 
giants. M dwarfs are identified by filled circles. Contrast values for 
the Sun are indicated by circled points and represent the umbral 
(higher) and penumbral (lower) temperature contrasts. Only the 
umbral contrast in included in the fits. Blue dashed line is linear 
fit to all points. Solid purple line is linear fit to dwarf stars only. 
Red dot-dash line is polynomial fit to dwarf stars only. Linear fit 
lines have limited use at predicting temperature contrasts at low 
photospheric temperatures due to the fact that they reach zero 
in the T = 2400 - 2600 K range. 

We developed the SPOTSS (Simulated Patterns Of 
Temperatures on Stellar Surfaces) code to simulate stellar 
spot patterns (see iKorhonen et al .1120141 . hereafter Paper I, 
for more details). Using SPOTSS we can generate either 
random spots, or spots constrained to a certain active lati¬ 
tude and/or longitude region. The spots created by SPOTSS 
have a umbral region which has a temperature defined as the 
spot temperature, and a penumbral region, with a temper¬ 
ature of halfway between the spot temperature and photo¬ 
spheric temperature. Spots are generated with a lognormal 
size dis tributio n (the s ize di stribution observed for Sunspots, 
see e.g. iBoedan et alii 19881) . combined with a power law to 
slightly increase th e number of small spots fo r M dwarfs due 
to observations bv lJackson fe Jeffriesl (|2012l l which indicate 
M dwarfs have numerous small spots covering the surface. 

In all our simulations we used a stellar rotation period 
of P = 7.7 days. Paper I examines the effects of vsin i on 
the measured jitter. We are interested here in exploring the 
effect of spot configuration on jitter and thus we keep a 
constant usin i for all our “observations.” A 7.7-day period 
is a reasonable choice for M dwarfs, and was chosen because 
it is in between the very active and the non-active cases. 
Since the choice of activity level (filling factor) and spot 
configuration have a much higher effect on the jitter than 
the usini, using P=7.7 d for all our measurements does not 
significantly bias our results. 

3.1 RV Jitter from spots 

Radial velocity measurements were generated f rom our 
simu l ated stellar surfaces u sing DIRECT7 (Pis kunov et al.l 
Il990l : lHackman et al j|2()()l| j and DEEMA (Detection of Ex¬ 
oplanets under the Effect of Magnetic Activity) codes. For 
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Table 1. Giant star spot temperature contrasts from literature. Where available, photometric 
variation (AV) and spot filling factor (ff) are also included. The method used to obtain the spot 
temperatures is indicated as: molecular band head modeling (MB), light curve monitoring (LC), 
line depth ratios (LDR), doppler imaging (DI), or eclipse mapping (EM). 


Star 

SpT 

Tphot 

Tspot 

AT 

AV 

ff (%) 

method 

reference 

HD199178 

G5III 

5350 

3800 

1550T200 


29 

MB 

1 * 

El Eri 

G5IV 

5600 

3700 

1900T150 


36 

MB 

1 * 

A And 

G8III 

4750 

3650 

1100T150 


23 

MB 

1 * 


G8IV 

4780 

3730 

1050T100 


17 

LC 

4 

HK Lac 

K0 III 

4790 

3840 

950T200 

0.14 


LC 

2 * 


K0 III 

4765 

3788 

977±10 

0.37 

34 

LDR 

3 * 


K0 III 

4600 

3520 

1080T100 


8 

LC 

4 


K0 III 

4765 

3565 

1200T100 

0.74 


LC 

5 


K0 III 

4765 

3998 

767T119 

0.07 

13 

LDR 

15 


K0 III 

4765 

4050 

715±91 

0.07 

14 

LDR 

15 

XX Tri 

K0 III 

4820 

3500 

1300 

0.63 

20 

DI 

6* 


K0 III 

4750 

3420 

1320T120 


35 

MB 

7* 

HU Vir 

K0IV 

5000 

3440 

1560T100 


44 

MB 

1* 

UX Ari 

K0IV 

5000 

3500 

1500 


48 

MB 

8* 


K0IV 

4780 

3360 

1420T100 


8 

LC 

4 

LX Per 

K0IV 

4780 

4050 

730T100 


5 

LC 

4 

AR Lac 

K0IV 

4700 

3500 

1200 

0.04 


LC 

9* 

a Gem 

K1III 

4600 

3850 

750T100 


33 

MB 

1* 


Kim 

4440 

3870 

570T100 


8 

LC 

4 

DM UMa 

Kim 

4600 

3570 

1030T100 


35 

MB 

1* 


Kim 

4500 

3450 

1010T120 


42 

MB 

7* 

SZ Psc 

K1IV 

4700 

3500 

1200T400 



LC 

10* 

HR1099 

K1IV 

4700 

3500 

1200 


40 

MB 

8* 


K1IV 

4700 

3500 

1200 

0.09 


LC 

9* 

HR7275 

K1IV 

4600 

3400 

1200T100 


9 

LC 

4 


K2III 

4600 

3500 

1100T150 


27 

MB 

1* 

IM Peg 

K2II 

4450 

3450 

1000 


20 

DI 

12* 


KIV 

4400 

3270 

1130T35 

0.38 

15 

LC 

13 


K2II 

4666 

4218 

448T126 

0.17 

11 

LDR 

15 


K2II 

4666 

4269 

397T140 

0.17 

12 

LDR 

15 


K2III 

4666 

3943 

723±10 

0.20 

32 

LDR 

3* 


K2III 

4440 

3520 

920T100 


7 

LC 

4 

II Peg 

K2IV 

4600 

3400 

1200T100 

0.42 

37 

LC 

2* 


K2IV 

4750 

3530 

1220T100 

0.16 

50 

MB 

1* 


K2IV 

4750 

3530 

1220T100 

0.23 

43 

MB 

1* 


K2IV 

4500 

3300 

1200 

0.25 


LC 

9* 


K2IV 

4600 

3500 

1100 

0.25 

20 

DI 

11* 


K2IV 

4350 

3620 

730T100 



LC 

4 

IN Vir 

K2IV 

4600 

3350 

1250T170 


40 

MB 

7* 

VY Ari 

K3IV 

4916 

4007 

909±10 

0.41 

41 

LDR 

3* 


K3IV 

4600 

3400 

1200T50 

0.28 

15 

LC 

14* 


K3IV 

4600 

3400 

1200T50 

0.20 

12 

LC 

14* 


K3IV 

4600 

3400 

1200T50 

0.10 

15 

LC 

14* 


K3IV 

4916 

4080 

836T118 

0.15 

15 

LDR 

15 


K3IV 

4916 

4129 

787±88 

0.15 

16 

LDR 

15 


References: 1, |()Nca 1 eta]] 1 1 : 2. IVogd ll 198 ll) : 3, [Catalano_et_al 
1 19851'! : ’ 

20011 ) 


01ali^t^^^9M)j6, Strassmeieij l[l99t ); 7~rO^Neal"Tl""aI 


Rodono etal.l i 19861) : 10 


Eaton & Hall 


fl97ai ~ 


11 , 


20031 : 

(200J); 


Berdvugina et aLTI^OOol i: 13. IPadmakar fe*^ande^lll999lh 14, l&t 

Frascaetlin~l2005ll : 


Berd^ugina i et al 


rassmeier 


fc Bopd 



* indicates that these values are also included in the iBerdvuginal ll2005ll review. 


M dwar f photospheric tempera tures we used the MARCS 
models dGustafsson et aDl2008l l. with solar metallicity, and 
logg = 5, a temperature grid of 2500 - 3900 K, with a 200 K 
step size, and 17 limb angles. DIRECT7 calculates synthetic 
spectra by integrating spectral line profiles over the entire 
visible stellar surface at each observational phase, using an 


evenly spaced wavelength grid. The vsin i and inclination 
values are taken into account to calculate which part of 
the stellar surface is visible and the rotational broadening 
of the spectral lines. The number of observational phases 
is evenly distributed over the length of the observing run. 
DEEM A obtains radial velocity (jitter) measurements by 
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Table 2. Dwarf star spot temperature contrasts from literature. Where available, photometric 
variation (AV) and spot filling factor (ff) are also included. 


Star 

SpT 

Tphot 

Tspot 

AT 

AV 

ff (%) 

method 

reference 

Sun 

G2V 

5870 

4200 

1670 


0.03 - 0.30 


25 


G2V 

5870 

5000 

870 




25 

EK Dra 

G2V 

5930 

5400 

500 


6 

LC 

16* 


G2V 

5850 

4800 

1050 

0.08 

11 

DI 

17* 


G2V 

5830 

3800 

2030 


40 

MB 

7* 

HD307938 

G2V 

5800 

3900 

1900 

0.06 

13 

DI+LC 

18* 

AB Dor 

K0V 

5200 

4000 

1200 

0.05 

12 

LC 

19* 


K0V 

5200 

3700 

1500 

0.12 

5 

LC 

19* 

LQ Hya 

K2V 

5175 

3650 

1525 


45 

MB 

8 * 

V833 Tau 

K4V 

4500 

3175 

1325 


45 

MB 

8 * 

EQ Vir 

K5e V 

4380 

3350 

1030T120 


45 

MB 

7* 

BY Dra 

MOV 

4100 

3500 

600T450 

0.08 

34 

LC 

2 * 


MOV 

4100 

3700 

400 


60 

LC 

20 * 


MOV 

4100 

3500 

600 

0.10 


LC 

9* 

YY Gem 

MOV 

3820 

3400 

400 



LC 

21 * 

KOI 1152 

Ml V 

4069 

3869 

200 



EM 

26 


Ml V 

4069 

4030 

39 



EM 

29 

HHJ 409 

MV 

3560 

3300 

260 

0.08 

13 

LC 

27 

AU Mic 

M2e V 

3500 

2650 

850 

0.10 


LC 

9* 


M2e V 

3500 

3000 

500 

0.32 

10 

LC 

22 * 

EV Lac 

M4e V 

3300 

3080 

220 

0.10 

7 

LC 

23* 

BPL129 

M4V 

3200 

2400 

800T200 


4 

LC 

24 

LHS 6351 

MV 

3150 

2900 

250 


1 

LC 

28 


References: 2, JVog^_|l98lj^; TjCTNeal et L jyjj[200^_8 ; _[O’Neal et aI]j[200l|V_9 : _ [Rodonoet al.1 
ll986|); 16jPorrea_^-GuiaaiJ ll994|); 17.IStrass meier fc Ricj (1199811; 18.lMarsden eta.l .1 j‘20Q. r J) : 19, 
|Amado_et_aL|j[2000h^0jChug^tao3 lll97(|kJ(l jTorrej^^4ibad ll^OO^kjgjTorres fc^terraT Mellol 
dl973h : 23. lAbranin et al.l lll998h ; 24~ [Scholz_et_alJ l|2005f); 25, [Berdmginal l2005|); 26. Sa nchis 
Ojeda , priva te communication; 27. lTerndrup"et al.l d 1 999l i : 28. ISavanovfc^ndtrienkol ll2012l i: 29, 
lYargal (2014i i 

* indicates that these values are also included in the iBerdvuginal (12005 1 review. 


cross-correlating the calculated line profiles with the line 
profile from the spectrum generated at the first observa¬ 
tional phase. See Paper I for a more detailed description of 
the codes and process. 

Figure [3] shows the RV jitter from calculations us¬ 
ing seven 30 A wavelength intervals (i.e. 4610-4640 A, 5510- 
5540 A, etc) ranging from 3710 A to 9110 A. Three differ¬ 
ent photosphere-to-spot temperature contrasts were investi¬ 
gated: a low-contrast case of AT = 100K, a medium contrast 
of AT = 500K, and a high contrast of AT = 900K. Our re¬ 
sults confirmed previous studies, which showed that jitter 
decreases with increasing wavelength; redder wavelengths 
have a lower spot-induced jitter. At high contrast, we saw 
a slight incre ase g oing from 8210 A to 9110 A, similar to 
IBarnes et al.l (|2011T ). Although it is tempting to conclude 
from our analysis that the best wavelength for planet hunt¬ 
ing around M dwarfs would be in the IR, it is also important 
to chose a bandpass with a large amount of spectral infor¬ 
mation for the cross-correlation to yield the most precise 
results. 

We also tested the RV jitter as a function of stellar pho- 
tospheric temperature (spectral type). We investigated the 
same low, medium and high AT values, as before, but we 
also added a 4th case where we used T spo t = 0.9 T e g for the 
spot temperature (e.g., for T e g = 3800 K, T apot = 3420 K) 
in order to keep the contrast ratio constant as we changed 


spectral type (see Figure [4J, since the photospheric temper¬ 
ature was changing and thus a constant AT would not result 
in a constant T e ff/T apo t ratio. We found that the RV jitter 
would first decrease, then reach a minimum, and then in¬ 
crease again with decreasing photospheric temperature (i.e. 
increasing spectral type). The approximate spectral type of 
this minimum was dependent on contrast ratio. At low con¬ 
trast, the effect was less extreme. It is unclear whether this 
effect is real or could be caused by possible problems with 
creating synthetic spectra at lower temperatures. 

With the exception of the Sun, we cannot directly ob¬ 
serve stellar spot filling factors. Therefore, we used DI¬ 
RECT? to investigate the photometric behavior calculated 
from the input spot configuration maps and how the vari¬ 
ation in magnitude relates to the filling factor, and the re¬ 
sulting RV jitter. Figure [5] shows the variation in AV mag as 
a function of spot filling factor, from 0.3% to 99.7% filling 
factor. The photometric variation increases with increasing 
spot filling factor to a peak around a filling factor of 50%, 
and then decreases again. This is expected since high filling 
factors of dark spots would simply appear as bright spots on 
a dark surface. The slope of this curve, however, is dependent 
on the size distribution of the spots, since a few large spots 
would make a larger difference in the photometric variation 
than many small, evenly-spaced spots. 

We created three different sets of spot configura- 
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Figure 3. Left: RV jitter generated by a single, 3-degree spot at five different observational wavelength ranges: 4610-4640 A, 
5510-5540 A, 6410-6440 A, 7310-7340 A, and 8210-8240 A. Right: Semi-amplitude of R.V jitter (normalized to the jitter measured 
at 3700 A) as a function of wavelength, at thee different photosphere-to-spot temperature contrasts: AT = 100 (triangles), 
AT = 500 (squares), and AT = 900 (circles). RV jitter is higher at all wavelengths for higher AT. Redder wavelengths reduce 
RV noise from spots, since the contrast between the photosphere and the spot is lower at these wavelengths, although this effect 
is less dramatic at low contrast ratios. 





AT(K) 


330 320 310 300 290 



Figure 4. Left: Semi-amplitude of RV jitter generated by a single, 10-degree spot at a range of photospheric temperatures from 
T e ff = 3800 K to 2800 K, corresponding to approximate spectral types M0 - M8. Three curves representing low (AT = 100, 
triangles), medium (AT = 500, squares), and high (AT = 900, circles) contrasts are shown. Jitter is normalized to the jitter 
measured at T e ff=3800K, AT = 900. Right: Semi-amplitude of RV jitter as a function of stellar temperature for a constant 
photosphere-to-spot temperature ratio given by T spo t/T e ff = 0.9 (e.g., for T e ff = 3800 K, T spo t = 3420 K). Jitter is normalized to 
the jitter measured at T e ff=3800K. 


tions based on G ranzer et al . I (iooO) an d I Vida et al.l I 2 OO 9 I) : 
Granzer , with random spots only emerging at an active lat¬ 
itude range, Vidal, with small random spots clustered in 
a distinct active region of both longitude and latitude, and 
Vida2, with two large spots separated by 150° in longitude, 
but varying in latitude (see Figure [6j). We also created a 
set of random spot distributions with a range of filling fac¬ 
tors. We compare the resulting jitter from each configuration 
in Figure [7] (left). The level of jitter, on average, increases 
with spot filling factor, however at a given filling factor the 
amount of jitter generated by a given spot map is highly 
dependent on the configuration of spots. 

Figure [7] (right) shows the RV jitter as a function of 
photometric variation of the star. Regardless of spot config¬ 


uration, the jitter increases linearly with photometric vari¬ 
ation. A line fit to all the data points is overplotted to il¬ 
lustrate this trend. Thus, stellar photometric variation is a 
much more useful predictor of RV jitter than spot filling fac¬ 
tor. This conclusion is suppo rted by observational evidence 
from, e.g. lLanza et al.l boill) . who used high-precision op¬ 
tical photometry by the MOST satellite to map the longi¬ 
tudinal distribution of active regions in late-t ype stars a nd 
to pre dict the RV jitter caused by the spots. Aigrain et alj 
J2012I) presented a method that predicted activity-induced 
jitter with high-precision time series photometry from a sim¬ 
ple spot model, and tested the method using MOST and 
SOPHIE observations of the planet host HD189733. Most 
recently, ICegla et al.l (120141 ') used both Kepler light curves 
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Figure 5. Observed V photometric variation of star with spot 
filling factor. Grey diamonds represent each individual measure¬ 
ment of jitter from a unique stellar surface consisting of randomly 
distributed spots with a given filling factor. Black diamonds cor¬ 
respond to the mean value of the measured jitter at each fill¬ 
ing factor, and error bars represent standard deviation calculated 
from all the measurements at each filling factor. Jitter increases 
to a maximum which is reached at approximately 50% spot cov¬ 
erage, and then decreases again as the spots dominate the stellar 
surface. 


Figure 8. Simulated radial velocity from planet + noise + jitter. 
20 “observations” (crosses) have been evenly distributed over an 
“observing run” with a length of 142 days. The planetary contri¬ 
bution to RV signal is overplotted (dashed line), and generated 
from a simulated system of a 0.5Mq star with a 1 Mq planet 
on a 58.7 day orbit. A single 3-degree spot on the stellar surface 
contributed to the jitter. The stellar rotation period is 7.7 days. 
We can see that the RV jitter from the spot is significantly higher 
than the “true” RV from the planet in many observations. 





Figure 6. Spot configurations, (a): Rand, random spots dis¬ 
tributed over the entire stellar surface, (b): Vidal, random spots 
at active longitude and lati tude, (c): Granzer, ra ndom spots at an 
active latitude range, as in lGranzer et all ll200d l. (d): Vida2, two 
large spots separated by 150° in longitude, similar to I Vida et al.1 

teoosl) . 


and GALEX data to examine the relation between RV jitter 
and photometric variation, concluding that for magnetically 
quiet stars the correlation was strong enough that photomet¬ 
ric measurements can be used to directly estimate jitter. 


3.2 RV from simulated planet 

We generated RV curves resulting from orbiting planets us¬ 
ing Kepler’s Third Law as is described in Paper I. Figure [8] 
shows the sinusoidal curve resulting from an orbiting planet 
(dashed line), and the RV measurements that would result 
from this curve plus the added jitter of a 3-degree spot on 
the stellar surface (crosses). We investigated some “detection 
limits” for orbiting planets by using a Lomb-Scargle Peri- 
odogram and attempting to recover the period of the planet 
from our radial velocity curves. These are not robust limits 
on planetary detection (see the following section) but simply 
a quick mechanism for analyzing the RV curves generated 
from the planet + spots model and how the noise introduced 
by the spots can obstruct planet detection. For each spot 
configuration + simulated orbiting planet, we investigated a 
number of different, evenly spaced sets of “measurements” 
to obtain an RV curve. We then created a Lomb-Scargle 
Periodogram of each RV curve and identified the planetary 
period in the resulting peaks. If there was a peak corre¬ 
sponding to the planetary period, we calculated the False 
Alarm Probability (FAP) of that period based on its Scar- 
gle Power. We calculated the FAP by randomly rearranging 
the RV data points N(N — 1) times (where N is a the num¬ 
ber of data points) and then plotting the randomized RV 
array against the time array, and taking a LS-Periodogram 
of the resulting curve. We then found any peaks in the peri¬ 
odogram, and their Scargle power, and compared these val¬ 
ues to the Scargle power of the original period detection. The 
FAP is the fraction of trials in which we found peaks with 
Scargle power higher than that of the peak corresponding 
to the actual period in the original, non-randomized data. 
This method of calculating the FAP is useful because it in¬ 
corporates the uncertainties already in the data, and does 
not require a gaussian noise distribution. 

We consider a FAP of ^ 1% to be a “detection” and 
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Figure 7. Left: Jitter from spots as a function of filling factor for the four different spot configurations shown in Fig urelfijr andom 
spots (green diamonds), spots at an active latitude (orange circles), spots in a distinct active region following IVida et al.l 120091 ') 
(light blue, downward triangles, Vidal), and two distinct spots separated by 150 degrees in longitude, at differing latitudes (dark 
blue, upward triangles, Vida2). Right: Jitter from the same spot configurations plotted as a function of V-band photometric 
variation of star. Whereas the jitter as a function of filling factor is highly dependent on spot configuration, when plotted as a 
function of photometric variation the resulting jitter shows almost no dependence on spot configuration but instead follows a 
linear relationship, which is shown by the overplotted dashed line (a linear fit to all data points). Vidal configurations show a bit 
more deviation from the fit line and also often resulted in significantly higher jitter for the same filling factors compared with the 
other three distributions. 


anything higher than this to be a non detection. We in¬ 
vestigated three contrast ratios, low, AT = 100, medium, 
AT = 500, and high, AT = 900, over a range of fill¬ 
ing factors, 0.1% (Solar) 0.3% (Solar maximum), 1%, 3%, 
10%, 30%, 50%, and 80%, with a “random” spot distribu¬ 
tion (see Figure [S] a). We simulated four different planets 
around each star with masses M = 1 A/® (Earth), 5 M® & 
10M® (Super-Earths), and 20M® (Neptune). We also in¬ 
vestigated two photospheric temperatures, T e ff = 3150 and 
T e ff = 3650, corresponding to stellar masses M* = 0.15 M® 
and M* = 0.50 M®, respectively. For each stellar temper¬ 
ature, we used a semi-major axis of orbit for the orbiting 
planets wi thin the classically define d conservative HZ for the 
star, from [Kopparapu et alj J2013I . see Section 14.31) . These 
orbits correspond to periods of 21.3 days and 84.7 days for 
T e ff = 3150 and T e fr = 3650, respectively. Unless otherwise 
specified, we used a default stellar rotation period of 7.7 
days. 

Figure [9] shows the FAP as a function of the number of 
measurements (“observations”) taken over a 112-day observ¬ 
ing run for the Solar filling factor cases of the T e // = 3650 
star with medium contrast ( T apo t = 3150, AT = 500). The 
solid green line represents the M = 1 M® planet and is not 
detectable even with 100 observations. The other planets are 
all detectable at these filling factors. 

“Intermediate” filling factors, 1% - 30%, are shown in 
Figure [IU1 As the filling factor increases, more planets evade 
detectability to a higher number of observations. At 1% fill¬ 
ing factor, the 1 M® and 5 M® are both undetectable with 
20 observations, but the 5 M® planet reaches a sufficiently 
low FAP of detectability at around 30 observations. By 30% 
filling factor, all of the planets are initially undetectable, 
although the 20 A/® planet reaches detectability with just 
a slight increase in number of observations, and the 10 M® 
planet is detectable in 70 observations. However, both the 


1M® and 5 M® planet remain undetectable even up to 100 
observations. 

“High” filling factors, 50% and 80%, are shown in Fig¬ 
ure [lT] Detectability is the lowest for all planets in the 50% 
filling factor case, and actually improves slightly at 80% fill¬ 
ing factor. The 80% filling factor case more closely resembles 
the 30% filling factor case shown in Figure flOl This makes 
sense since a star with an 80% coverage of dark spots is es¬ 
sentially the same as a dark star with 20% coverage of bright 
spots (see Figure [5j. 

The effect of temperature contrast on detectability is 
shown in Figures [12] and [T3] Figure \12\ illustrates the T e fr = 
3650 (M* = 0.50 Mq) case. The Solar filling factor (0.1%) 
case is shown on the left with low (a), medium (b), and high 

(c) contrast ratios. For the low contrast case, all planets are 
detectable. The 1 M® planet is not detectable for medium 
or high contrast at this filling factor. The plot on the right 
in Figure [12] shows the 30% filling factor case, again at low 

(d) , medium (e), and high (f) contrasts. Increasing contrast 
inhibits detectability of all planets, requiring more observa¬ 
tions to yield a detection. At high contrast, only the 20 A/® 
planet achieves detectability in this observational scenario. 
Figure [13] shows the same cases, but for the T e g — 3150 
(M = 0.15 Mq) star. A similar trend is observed, although 
detectability at a given planet mass + contrast + filling fac¬ 
tor case is significantly better around this star, due to the 
combination of smaller stellar mass and closer-in HZ orbit. 
Since we used the same observing period of 112 days for 
both stars, the data covers more orbits for the smaller star, 
since the closer-in HZ means a shorter orbital period. When 
the data are phase folded, more of the activity noise can 
be averaged out because at a given planetary phase there 
are more data points. We repeated the experiment using an 
observing period of 445 days for the 3650 K M dwarf (with 
a HZ orbit of 84.7 days), which corresponds to the same 
planetary orbital phase coverage as the 112-day observing 
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period of the 3150 K star (HZ orbit of 21.3 days). In these 
observations the detection limits were more similar to the 
3150 K case due to the similar phase coverage, although still 
not identical because the stellar RV generated by a planet 
of a given mass is different in the two systems due to the 
different stellar masses and HZ orbits corresponding to the 
different stellar temperatures. 

At a given filling factor, the spot contrast makes a huge 
difference in the number of observations required to detect 
these HZ planets. Thus, better constraints on spot tempera¬ 
tures would help to better plan observing strategies when 
searching for HZ planets, especially around active stars. 
Even at low filling factors, the difference in spot tempera¬ 
ture could mean the difference between finding a HZ planet 
and missing it. In Figure (low contrast case) the 1M® 
planet would be detected even with a low number of obser¬ 
vations. However, when the contrast is increased (112 b ). the 
FAP for this planet is still just above the detection thresh¬ 
old, even at 100 observations. When the contrasts is highest 
(HO, the FAP is an order of magnitude over the detection 
limit with 100 observations, and would require significantly 
more observations to achieve detectability. We extended the 
observations of the 1 M® planet (orbiting a 3650 K star with 
a 0.1% spot filling factor), and found that even at 500 ob¬ 
servations the FAP still did not drop below the detection 
threshold of 0.01. 

We also tested the effect of orbital eccentricity by inves¬ 
tigating a 5 M® planet with four orbital eccentricities rang¬ 
ing from 0.0 to 0.8. Figure [lo] illustrates a high (30%) and 
low (0.3%) filling factor case, and the same trend is seen 
in both: higher eccentricity causes a lower FAP (though the 
two highest eccentricity cases, 0.5 and 0.8, resulted in almost 
the exact same FAP for all measurements). Thus, a planet 
with a more circular orbit would be more difficult to detect 
than a planet with higher eccentricity if all other stellar and 
planet parameters were equal, and if the signals are well- 
sampled in phase. This is due to the fact that a planet in a 
highly eccentric orbit will create an RV signal with a higher 
semi-amplitude, as seen in Figure 1141 However, if the sig¬ 
nal is not well-sampled, the signature of a high-eccentricity 
planet could be missed completely, because the RV curve is 
mostly relatively flat, with short-lived, high-amplitude peaks 
(or troughs). 


4 UNCERTAINTY IN DERIVED PLANETARY 
SYSTEM PARAMETERS 

We hesitate to declare hard “limits” on exoplanet detection 
due to the ever-increasing precision with which RV curves 
can be measured and the many techniques that are being 
constantly developed to improve RV data and remove spuri¬ 
ous periodic signals in order to uncover planetary signatures, 
such as Fourier component analysis v ia pre-whitening, an d 
local trend filtering (see, for example, iBoisse et al.l (1201 lblj . 
who corrected for activity noise by fitting sinusoids at the 
stellar rotation period and its two-first harmonics and were 
able to remov e about 90% of the RV jitter amplitude; and 
lHatzesI d2013h . who demonstrated that it is possible to ex¬ 
tract the RV signal of a short-period Earth-mass planet de¬ 
spite stellar activity noise, if high-cadence observations are 
available). Developing technology is allowing spectra to be 



Figure 14. Stellar radial velocity curves generated by orbiting 
planets with identical parameters but different orbital eccentric¬ 
ities: 0.0 (solid line) and 0.9 (dashed line). Y-axis is in arbitrary 
velocity units where the high eccentricity case has been normal¬ 
ized to 1. X-axis is also normalized to one in arbitrary time units. 



Figure 15. LS FAP of planet detection as a function of number 
of measurements for a 5 M® planet with four different orbital 
eccentricities: 0.0 (red solid line), 0.3 (dark orange, dashed line), 
0.5 (orange dash-dot line), 0.8 (yellow dotted line). FAP detected 
threshold is indicated by dashed grey line. All orbits have the 
same period, 84.7 days. The stellar temperature is 3650 K and the 
spot temperature is 3150 K (AT = 500 A). Spot filling factors of 
30% (top) and 0.3% (bottom) are compared. The same trend is 
seen in both: high eccentricities have lower FAP. The 0.5 and 0.8 
eccentricities, however, are almost exactly the same. 


measured to such precision that the limit on derived radial 
velocities now stems from li mits due to properties of the 
star itself. iDumusaue et all (I2OI2I ) analyzed high-precision 
radial velocities for a Centauri B, and identified the follow¬ 
ing contributions to the radial velocity signal: instrumental 
noise, stellar oscillation modes, granulation at the surface of 
the star, rotational activity, long-term activity induced by 
a magnetic cycle, the orbital motion of the binary system, 
light contamination from a Centauri A, and imprecise stel¬ 
lar coordinates. They then attempted to model and remove 
each signal separately. They extracted the RV sign ature of 
a small, approximately Earth-mass planet. However. lHatzesI 
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Figure 9. Lomb Scargle (LS) false alarm probability of planet detection as a function of number of measurements made during 
a 112-day period. The two Solar filling factor cases: 0.1% (left) and 0.3% (right) are shown. Stellar temperature is T e ff = 3650 
and stellar mass is M* = O.15M0. Spot temperature contrast is AT = 500, for s spot temperature of T spo t = 3150. Planets 
are on an 84.7-day orbit which corresponds to an orbital radius within the HZ for this star. Four planet masses are included: 
M = 1 M® (solid green line), 5M® (dashed blue line), 10M® (dot-dashed purple line), and 20 M® (dotted red line). The M 
= 1 M® remains above the detectability (defined as a FAP of 1% or less) threshold (indicated by the dashed grey line) for all 
numbers of measurements. The 5 M® planet is also detectable in all instances, although at 20 observations the FAP just grazes 
the detectability limit in the 0.3% filling factor case (right). 10 M® and 20 M® planets have a FAP of zero consistently in both 
cases. Note that zero is plotted at a line just above the x-axis. 
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Figure 10. LS FAP of planet detection as a function of number of measurements for the same star and planet cases as in Figure 
[9] Four “intermediate” filling factor cases are illustrated: 1% (top, left), 3% (top, right), 10% (bottom, left) and 30% (bottom, 
right). Increasing the filling factor from 1% to 10% (left) requires an increase from 30 to 80 measurements for the 5M® planet to 
be detected. Grey dashed line indicates FAP detection threshold: planets are considered “detectable” where FAP falls below this 
line. The same trend is seen with the 10 M® planet (right) where an increase from 3% filling factor to 30% requires a jump from 
20 to 70 measurements to reach the detectability threshold. 

(12013ft reanalyzed the data using different methods to re¬ 
move the stellar activity signal from the RV data, and was 
unab le to recover a s ignific ant RV signal corresponding to 
the lDumusaue et al.l d2012ft planet. 

Thus, we do not attempt to set “limits” on radial ve¬ 
locity measurements based on this study alone. Rather, in 
this work we are interested in the potential uncertainty in 


planetary system parameters derived from RV observations 
that can be introduced via contributed stellar jitter. 

To derive plane tary parameter s from our synthetic RV 
curves, we used the lEastman et all d2013ft EXOFAST suite 
of IDL routines. EXOFAST was created to fit exoplanetary 
transits and radial velocity variations either simultaneously 
or separately, and characterize the parameter uncertainties 
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Figure 11. LS FAP of planet detection as a function of number of measurements for the same star and planet cases as in Figure 
[9] “High” filling factors of 50% (left) and 80% (right) are illustrated. The number of observations required to achieve detectability 
for all planets decreases slightly in the 80% filling factor case from the 50% case. The 80% case more closely resembles the 30% 
case from Figure [To] This trend was expected since with spot filling factors of above 50% the spots dominate the stellar surface 
and the level of spot-induced jitter actually decreases. 
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Figure 12. LS FAP of planet detection as a function of number of measurements for a T e g = 3650 star. Solar filling factor (0.1%) 
case is shown on the left at low (a) medium (b) and high (c) temperature contrasts, and 30% filling factor is shown on the right, 
also at low (d), medium (e), and high (f) temperature contrasts. 1M® planet (solid green line) is only detectable in low filling 
factor, low contrast case (a). Note that in the (a) panel of the left plot, FAP is 0 for all planet masses. 
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Figure 13. LS FAP of planet detection as a function of number of measurements for the same filling factor and temperature 
contrast cases as in Figure fl2l but for a T e ff = 3150 star. In panels (a), (b), and (d) FAP is 0 for all planets. 
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and covariances using a Differential Evolution Markov Chain 
Monte Carlo method. Since we are interested in RV jitter, 
we only fit the RV variations with EXOFAST. Simultaneous 
fitting of transits and RV variations will be investigated in 
a later work. 

Markov Chain Monte Carlo (MCMC) algorithms are 
commonly used to fit data points to a model with mul¬ 
tiple parameters. EXOFAST uses the Metropolis-Hastings 
MCMC algorithm to sample the probability that a given 
model M is correct, given the data, D: P(M\D), starting 
with a trial set of model parameters (see below for a discus¬ 
sion of the specific parameterization we used for our models) 
and evaluating y 2 with respect to the data. EXOFAST then 
calculates the ratio of the likelihood of a new set of pa¬ 
rameters relative to the initial set. Then a uniform random 
number is drawn from [0,1]. The new model is rejected if 
the random number is greater than the likelihood ratio, and 
the previous model is duplicated and added to the “chain.” 
If the random number is smaller than the likelihood ratio, 
the model is accepted, added to the chain, and the process is 
repeated, stepping to another region of parameter space. Be¬ 
cause we are assuming that the RV curve is generated from a 
star-planet system and can be described by Kepler’s Laws, 
we do not consider the likelihood of models with different 
parameterizations. 

EXOFAST improves on the MCMC method by imple¬ 
menting Differential Evolution MCMC (DE-MC). DE-MC 
is a more elegant approach to the somewhat “brute force” 
Monte Carlo method of parameter fitting. While the regu¬ 
lar MCMC method runs a single chain, the DE-MC method 
will run a number of chains equal to twice the number of 
free parameters in parallel and calculates the next step in 
the chain by taking the difference between the parameter 
values between two random chains. This approach will dra¬ 
matically decrease the number of links it takes for the chains 
to become “well-mixed” and to sufficiently sample param¬ 
eter space. Since it is not practical to use the same step 
size for many different parameters with different units, EX¬ 
OFAST adds the additional modification of creating self- 
adjusting step-sizes by varying each parameter individually 
until Ax 2 = 1 and then adding a uniform deviate equal to 
1/10 of that step size. 

The parameterization for RV fitting that EXOFAST 
uses is logP, log K, y/ecos w*, y/esinuj*, Tc, 7 , and 7 , where 
P is the period of the orbit, K is the RV semi-amplitude (in 
m/s), e is the eccentricity of the orbit (see section roii , w* 
is the argument of periastron of the orbit, Tc is the time 
of central transit (not used in RV-only fits), and 7 is the 
RV systemic velocity (zero point), in m/s. The fit will al¬ 
low us to obtain K, e, uj, P, Tp, Tc, as well as derived pa¬ 
rameters such as M t , R *, msini, and a. One alteration we 
made to the standard release of the EXOFAST code was to 
change t he estimate o f the s tellar mass/radius. EXOFAST 
uses the I Torres et al.l (120091 ') relation, which was not de¬ 
signed to work for cool stars with mass < ~ 0, 6 M r . We 
replaced the Torres relation with the [Baraffe et all ( 19981 ') 
models in order to corre ctly estimate the der ived parameters 
M*, /?*, msin i, and a. iBaraffe et al.l (1 19981 1 agrees reason- 
abl y well with Delfossg_et_aL ( 2000l l . and was chosen over 
the lDelfosse et a l.| (2000 |) for consis tency, because DEEM A 
implements the Baraffe et al.1 (119981 1 models in order to com¬ 
pute stellar radii (see Paper I). Since K, e, and P do not 


depend on the stellar mass/radius relation, they are inde¬ 
pendent of which model we choose. 

By default, EXOFAST handles uncertainties by scaling 
the uncertainties by a constant factor. A jitter term is not 
included (a systematic such as jitter that is independent of 
the signal should be added in quadrature), which is a com¬ 
mon practice with RV fits. Although we considered editing 
the code to include a jitter term, which might improve the 
quality of the fits, it seemed more useful to carry out analy¬ 
sis of the RV data in the usual way, rather than to account 
for a systematic error that we are aware of but observational 
planet hunters would likely not be. Thus, we analyzed our 
data in the standard EXOFAST procedure. 

EXOFAST allows the maximum and minimum period 
for the Lomb-Scargle periodogram to be set, in order to 
investigate periods that are slightly disfavored in the peri¬ 
odogram. We used this option to eliminate the rotational 
period of the star from the results, since this tends to domi¬ 
nate the signal. EXOFAST r eturns a large variet y of outputs 
from the fit parameters fsee [Eastman et aLll2013l l. The pa¬ 
rameters used in this work are: e, the eccentricity of the 
orbit; P, the period of the orbit, in days; and a, the semi¬ 
major axis of the orbit, in AU. 


4.1 Period 

We chose a 10 M© planet with a period of 21.2 days or¬ 
biting a 0.57 Mq M dwarf with a stellar rotation period of 
7.7 days. This period corresponds to an orbital radius that 
would put the planet wel l inside the classically-d efined hab¬ 
itable zone, according to dKopparapu et aT]l2013l . see Section 
\U\ for further discussion). We used five filling factors, 0.1% 
(Solar), 0.3% (approximatly the Solar spot filling factor at 
Solar maximum), 2%, 10%, and 30%. Tests with filling fac¬ 
tors above ~30% resulted in too many of the trials being 
unable to extract the planetary RV, and thus we were not 
able to make a meaningful statement about the uncertainties 
in the derived planetary parameters for these higher filling 
factors. Example maps of four filling factor cases are shown 
in Figure 1161 For each filling factor, we created 30 random 
spot maps for each of three different spot size distributions 
resulting in a total of 90 spot maps per filling factor. We 
then ran each spot map through DEEMA, which calculated 
the RV jitter contribution from the star and added to that 
the RV created by the planet’s orbit at each observational 
phase. The result was a curve of radial velocity (including 
the planetary contribution plus jitter) as a function of time. 

We took the resulting RV curves and reduced each one 
with EXOFAST, in an attempt to recover the planetary 
system parameters. We used the maximum- and minimum- 
period input keywords to eliminate the stellar rotation pe¬ 
riod, which tends to be the strongest periodic signal in the 
data due to the highly-spotted stellar surfaces. At each fill¬ 
ing factor, we calculated the mean planetary period, the 
standard deviation, and the standard error of the mean. We 
also ran an unspotted surface through the same process for 
calibration, and this test case is labeled as filling factor of 
zero. There is no standard deviation or expectation value 
for this case since we only used one map. We found that 
higher spot filling factors result in higher uncertainty in the 
final measured planetary period (see Figure □3). As filling 
factor increases, the measured values for the period deviate 
































M dwarf activity as noise in exoplanet detection 15 



Figure 16. Examples of stellar surface with five different filling 
factors, (a): 0.1%, (b): 0.3% (approximately the active solar filling 
factor), (c): 2%, (d): 10%, (e): 30%. 

increasingly from that actual value, although the mean value 
remains close to the actual period until the 30% filling fac¬ 
tor case. Even with small filling factors, e.g., 0.3%, some of 
the period measurements are very different from the input 
period. 

As the spot filling factor increased, the percentage of 
tests where the planetary period was not recovered from the 
RV curve also increased. For the Solar filling factors, 0.1% 
and 0.3%, the planetary signature was extracted in all tests. 
For filling factors of 2%, 10%, and 30%, the percentages of 
non-detections were 30%, 50%, and 60%, respectively. Be¬ 
cause larger filling factors resulted in many RV curves where 
the planetary signature is undetectable above the noise from 
the RV jitter, we focused on filling factors of 30% and be¬ 
low, although we point out that some active M dwarfs are 
thought to have filling factors of 50% or higher. 

4.2 Eccentricity 

The choice of eccentricity to use into our planetary system 
models was rather arbitrary, and we chose an initial eccen¬ 
tricity of zero for two reasons: (1) a circular orbit is the 
most simple case; and (2) from our investigation of FAP 
we found that, with a signal well-sampled in phase, planets 
with higher orbital eccentricities are easier to detect in a LS 
periodogram. We chose to use a circular orbit which adds 
the most uncertainty to planet detection using our sampling 
methods. 

We used the same spot maps as in the investigation of 
period, but instead looked at the derived eccentricity of the 
planetary orbit from EXOFAST. The results from the eccen¬ 
tricity study are shown in Figure [18] Similar to Figure [T 7 ] 
we calculated the mean eccentricity, the standard deviation, 
and the standard error of the mean for all the measurements 
at each filling factor. When the filling factor increases, the 
measured values for the eccentricity deviate substantially 
from the actual value. For filling factors of 10% and 30% 
basically all the measurements give large eccentricities of 
> 0.5, and even a relatively small filling factor of 0.3% can 
yield eccentricity measurements as high as 0.5. We acknowl¬ 
edge that at the ’’zero filling factor” point (filled diamond), 
i.e. the unspotted star, a greater than zero eccentricity was 



Figure 18. Eccentricity derived from radial velocity curves us¬ 
ing EXOFAST, at four spot filling factors. Filled circles indicate 
each individual eccentricity value derived at corresponding spot 
filling factor. Dashed line indicates zero eccentricity, the eccen¬ 
tricity value used in generating the simulated RV due to an or¬ 
biting planet. Open diamonds correspond to the mean value of 
all the measured eccentricities at each filling factor. Red error 
bars show the standard deviation of the measured eccentricities. 
Shaded gray region indicates the standard error of the mean. 


also measured, likely due to the Lucy-Sweeny Bias (see be¬ 
low). Still, the large measured eccentricities seen with filling 
factors 0.3% and larger cannot be explained by the Lucy- 
Sweeney bias. It is interesting to note that many exoplan- 
ets show larger eccentricities than th e Solar system planets 
llButler et al]|2006l ; iKane et akll2012l ), and one possible ex¬ 
planation could be spot activity of the host star. As was 
shown in Paper I, a spot located at the ‘right’ place can in¬ 
troduce changes in the measured radial velocity curve that 
can be interpreted as larger eccentricity. 


4-2.1 Eccentricity bias 

There exists a long-understood observational bias against 
low eccentricities in binary systems due to the data re¬ 
duction methods of observations of such systems, called 
the Lucy-Sweeney Bias. Systems with a circular orbit 
will be found to have a small but nonzero eccentricity 
due to observational u ncertainties and analytical errors. 
iLucv fe Sweenevl (Il97ll ) described this bias in binary star 
systems and how it results from a lack of phase space at 
exactly e = 0, so any uncertainly at all will push the re¬ 
sult toward a positive value, since ther e is no such thing as 
a negative eccentricity. lEastman et alj (j2013i ) has dramati¬ 
cally reduced the Lucy-Sweeney bias by allowing a negative 
eccentricity to be inferred, although they still caution that 
it is up to the user to inspect the PDF to determine the 
significance of the calculated eccentricity, and note that or¬ 
bits with intrinsically small uncertainties will have an over¬ 
estimated uncertainty with an uncertainty of roughly 0.007. 
Since we use an eccentricity of zero, we expect that the ma¬ 
jority of our eccentricities would be slightly overestimated, 
even without the added RV jitter. 
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Figure 17. Left: Derived periods and uncertainties as a function of spot filling factor. Periods were derived from radial velocity 
data (spot jitter plus RV from a simulated planetary system) using EXOFAST. Filled circles indicate each individual period 
derived at a given spot filling factor. Dashed line indicates the period of the planet used to generate the RV cure. Filled black 
diamond corresponds to “zero filling factor” case, i.e., an unspotted star. Open diamonds correspond to the mean value of all the 
measured periods at each filling factor. Red error bars show the standard deviation of the measured periods. Shaded gray region 
indicates the standard error of the mean. Right: Fraction of nondetections of the orbiting planet for each filling factor. 


4.3 Orbital radius and habitability 

The so-called “Habitable Zone” around a star is defined as 
the distance from the star where an orbiting planet could 
support liquid water on its surface. It is also sometimes 
referred to as the “Goldilocks Zone” since it is “not too 
hot, and not too cold, but just right.” Though this is a 
very Earth-centric definition of habitability, it is assumed 
(and mostly agreed upon) that some form of liquid is nec¬ 
essary for the evolutio n of lif e , and liquid water seems to 
be the most prevalent. ISeaeeij (I201fj j recently redefined the 
HZ taking into account to the possibility of thick planetary 
atmospheres trapping heat generated from active planetary 
interiors, thus extending the HZ farther from the star than 
was previously expected (and even allowing for the possibil¬ 
ity of habitable so-called “rogue planets”—planets that are 
not in a stellar orbit. However, for the purposes of this paper 
we will ado pt the more classic defin ition of the HZ, and use 
values from lKopparapu et al . ] rt2013h . which allow for both a 
“conservative habitable zone” and an “optimistic habitable 
zone”. 

The optimistic HZ is bounded by the “recent Venus 
limit” on the inner edge and the “early Mars limit” on the 
outer edge. These limits are empirical, based on the assump¬ 
tion that Ve nus has not hosted liqu id water for a least 1 
billion years dSolomon fe Headlll99ll ) and the evidence that 
Mars had liquid water on its surface about 3.8 billion years 
ago. Calculating the solar luminosity at these times com¬ 
pared to the present value, a HZ limit can be derived. More 
conservative HZ limits rely on the “water loss” (inner) and 
“maximum greenhouse” (outer) limits, which are simply, ac¬ 
cording to models, the highest stellar irradiance where all 
liquid water will not be lost, and the lowest stellar irradi¬ 
ance where the greenhouse effect could still keep a planet’s 
surface sufficiently warm to have liquid water, respectively. 

Whether or not a planet lies in the HZ depends on the 
stellar irradiance the planet is receiving, which is highly and 
sensitively depende nt on the star’s T e tf. To calculate an in¬ 
put value using the iKopparapu et ahl d2013l ) HZ equations 


(also available as an on line HZ calculator ) , we u sed the M 
dwarf luminosity from ICasagrande et ahl (l2008h . choosing 
the M dwarf HIP897 that has an effective temperature of 
3786 ± 126 K and a bolometric magnitude of 9.628 mag. Us¬ 
ing this valu e we calculated L t /Lp, = 0.011. Feeding this 
value into the lKopparapu et al.1 d2013ll script resulted in the 
following HZ limits: conservative, 0.1128 - 0.2102 AU, op¬ 
timistic, 0.0853 - 0.2197 AU. (See Figure [T9]) It would be 
naive to assume that planetary habitability could be com¬ 
pletely determined from total stellar irradiation alone, as 
there are many other factors that contribute to a planet’s 
habitability. Since M dwarfs have a very different spectral 
energy distribution from solar-type stars, planets orbiting M 
dwarfs will experience a distinct energy environment even 
at the same stellar irradiation as planets around a G dwarf 
dTarter et al.ll2007l ). The albedo feedback is dependent on 
the reflectivity of water/ice surfaces, which is wavelength- 
dependent and will be significantly differ ent for_ M dwarfs 
than f or the Sun. Also, the calculations of lKopparaou et al.l 
d2013l) do not account for the effects of water clouds, and 
do not apply to planets that are highly unlike Earth, either 
in atmospheric composition or surface gravity. These, along 
with many other factors compel us to present our “habitable 
zone” results with some caution. 

We used our simulated planet of P = 21.2 days around 
a 0.57 Mq M dwarf with eccentricity of zero, which results 
in an orbital radius of a = 0.12AU. This is well inside the 
conservative HZ for this star. Figure [T9] shows the measured 
a for the planet’s orbit as a function of spot coverage (filling 
factor), to demonstrate the uncertainty jitter can contribute 
to these measurements. It can be seen that there are cases 
where the measured value of a falls outside the habitable 
zone, although the mean and a for most filling factors still 
lie within this zone. Even at a filling factor as low as 2%, 
some of the cases resulted in planetary orbits outside the 
HZ. With increasing filling factor, the amount of scatter 
in the derived planetary orbits also increases. We did not 
include filling factors above 30% because at that level of 
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Figure 19. Derived semi-major axis (a) of orbit (filled circles) 
at each of four filling factors, plus “zero filling factor” which in¬ 
dicates an unspotted surface (filled diamond). Open diamonds 
indicate mean values of all measured a values at corresponding 
filling factor. The standard deviation of the measured values is 
displayed as red error bars, and the standard error of the mean 
is ind icated by a shaded g ray region. Habitable zone estimates 
from lKopparapu et alj J201 j l are included. The conservative HZ 
(indicated by green dotted lines) is defined by the ’’Moist Green¬ 
house Limit” and the ’’Maximum Greenhouse Limit”, and the 
optimistic HZ (indicated by blue dot-dash lines) is bracketed by 
the ” Recent Venus Limit” and the ” Early Mars Limit”. 

spot coverage too many tests resulted an undetectable plan¬ 
etary signal. However, for the cases where the planet was 
still detectable the scatter from the true a value was gener¬ 
ally higher, indicating that this increasing trend continues. 

As mentioned previously, high filling factors resulted in 
high percentages of non-detections for the planetary signal. 
In the cases where the planet was detected, we calculated the 
probability that this HZ planet would be classified as out¬ 
side the conservative HZ, P(O). For the Solar filling factors, 
0.1% and 0.3%, P(O) = 0, meaning all the planet detections 
were classified as habitable. The scatter in the measured a 
values due to jitter, was minimal although still present. For 
the higher filling factors, 2% - 30%, P(O) = .23. Thus, in 
approximately 20% of cases the HZ planet was classified as 
outside the HZ. 

All these tests were performed with a planet of mass 
10 M®. Lower planetary masses would mean the jitter sig¬ 
nal would be higher in comparison to the true RV, and thus 
we would expect a higher percentage of non detections and 
more scatter in the derived parameters for a given filling 
factor. A more extensive exploration of parameter space in¬ 
cluding many different eccentricities, planetary masses, stel¬ 
lar inclinations, and spot filling factors is the subject of a 
future study. 


5 CONCLUSIONS 

In this work we have investigated the effects of activity- 
induced jitter on planetary system parameters derived 
from RV fitting. We modeled the RV jitter using realistic 
spot models, and improved on previous jitter studies by 
including active latitude and longitude ranges, evidence 


of which has been seen in various observations of spot 
distributions on low-mass stars. Since published values of 
photosphere-to-spot temperature contrasts for M dwarf 
stars are still highly uncertain, we chose to investigate 
low, medium, and high contrast spots spanning the range 
of detected measured temperature contrasts from the 
literature. We quantified the error introduced to derived 
planetary parameters as a result of RV jitter, focusing on 
planetary period, eccentricity, and semi-major axis of orbit. 
We reached the following conclusions: 

• Jitter increases with increasing spot contrast at M 
dwarf photospheric temperatures, which confirms previous 
studies. This effect is most prominent at ~5500 A where the 
jitter at low contrast (AT = 100) is approximately a tenth 
that of the high contrast (AT = 900) case. 

• RV jitter decreases when moving to longer wave¬ 
lengths, as has been seen in previous studies. At high con¬ 
trast, when moving from 3710 A to 8210 A the jitter is re¬ 
duced to ~20% of the jitter at 3710 A. This effect is less 
extreme at low contrast, but even at AT = 100 the jitter at 
8210 A decreases to half that at 3710 A 

• When the spot contrast ratio is kept at a constant 
T a p 0 t/T e ff = 0.90, the RV jitter is lowest at earliest spectral 
types (i.e. highest photospheric temperatures). 

• At a given filling factor, the RV jitter from differ¬ 
ent spot distributions varies significantly. However, RV jitter 
increases monotonically with increasing photometric varia¬ 
tion, regardless of the spot configuration. In our tests the 
relationship was linear. Thus it is more useful to use the 
spot-induced stellar photometric variation than spot filling 
factor when trying to predict the RV jitter. Spot filling fac¬ 
tor is not a good predictor of spot-induced jitter. 

• The false alarm probability of planet detection in¬ 
creases with increasing filling factor as well as with increas¬ 
ing contrast between spot and photosphere. We found that 
an Earth-sized planet orbiting in the HZ of a T e // = 3650 K 
star with a solar filling factor of 1% could not be detected 
except in the lowest contrast case. With a filling factor of 
30% the same star could easily mask planets of 1 and 5 
Earth-masses with medium or high contrast, and even at 
low contrast the 1 M® is never detectable even up to 100 
observations. In the HZ of a T e // = 3150 if star, the detec¬ 
tion ability was increased. The Earth-sized planet could be 
detected in a reasonable number of observations at all con¬ 
trasts for the 1 filling factor case, and for 30% filling factor 
the Earth-sized planet was detectable at low contrast. 

• Higher orbital eccentricities of the planet result in 
lower FAPs, meaning planets on a circular orbit are more 
difficult to detect in the presence of RV noise than those 
on a more eccentric orbit. However, the difference in FAP 
moving from an eccentricity of 0.5 to 0.8 is practically zero 
in all cases, while the difference between the FAP at 0.0 
eccentricity and at 0.3 eccentricity can be quite large (up to 
2 orders of magnitude), especially at high filling factors. We 
note that this will only consistently be true if the RV signals 
are well-sampled in phase. 

• For tests with spot filling factors up to 10%, the mean 
and expectation value of the derived planetary period is still 
close to the actual value. However, at a 30% filling factor 
both the mean and expectation value deviate significantly 
from the actual period, although the actual period still lies 
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within the standard deviation of all the measured values at 
this filling factor. At filling factors higher than ~30%, a sig¬ 
nificant number of the tests could not recover the planetary 
period at all. 

• Due to the Lucy-Sweeny bias, we were not expecting 
to recover the eccentricity of 0 in the original planetary or¬ 
bit. Indeed, even the “zero filling factor” star, or in other 
words, the unspotted surface, resulted in an eccentricity of 
slightly above zero. However, as the filling factor increases, 
the measured eccentricity radically increases as well, indi¬ 
cating that with highly spotted stars, one should be wary of 
a measured eccentricity value. 

• At high filling factors, it is possible for the derived 
planetary orbits to scatter outside of the conservative HZ, 
although the mean and expectation values measured at each 
different filling factor still remain well within these limits. 
At filling factors of 2% - 30% the probability of classifying 
the HZ planet as orbiting outside the conservative HZ is 
23%. (However, we caution that orbital radius alone is not 
sufficient to truly determine habitability, and that there 
are still many factors that are not taken into account when 
calculating the HZ distance.) 

We have shown that stellar RV jitter induced by ac¬ 
tivity can contribute an important source of noise to plan¬ 
etary measurement that must be accounted for both when 
searching for and when characterizing potentially habitable 
planets around active stars. Better observational constraints 
on M dwarf spot parameters-spot sizes, distributions, and 
temperatures-will help to further constrain the effects of 
RV jitter on planet detection and characterization. Stellar 
photosphere-to-spot temperature contrasts seem to be the 
most important factor in predicting RV jitter. Large tele¬ 
scopes with high-resolution spectrographs will permit more 
Doppler Images of M dwarf spots to be taken in the future, 
and photometric monitoring from missions such as Kepler 
and TESS should provide further opportunities to constrain 
spots using eclipse mapping. Predicting RV jitter based on 
photometry alone removes the need to know the actual spot 
configurations, which is difficult for M d warf stars d ue t o 
their low lum i nosities. Studies such as iLanza et ali (1201 ill . 
lAigrain et al.l d2012lf ICeela et al.l (120141 ') , and this work have 
demonstrated the strong correlation between photometric 
variation and radial velocity jitter, which makes such an ap¬ 
proach promising. 
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